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Abstract 

In this paper we present a comprehensive theoretical description of molec- 
ular spectral converters in the specific context of Luminescent Solar Con- 
centrators (LSCs). The theoretical model is an extension to a three-level 
system interacting with a solar radiation bath of the standard quantum the- 
ory of atomic radiative processes. We derive the equilibrium equations of 
the conversion process and provide specific examples of application of this 
principle to the development of solar concentration devices. 



1. Introduction 

Luminescent solar concentrators (LSCs) have been first introduced in the 
late seventies as one of the simplest methods of concentrating sunlight in 
thin polymeric slabs [Tj and have been lately reconsidered in the light of 
the recent advances in material science and nano-optics as a cheap alterna- 
tive to standard photovoltaic modules [21 [3l E] . The concentration process 
can be understood through the principles of classic geometric optics. The 
polymeric slabs are doped with active molecules, which absorb a portion of 
the incident solar radiation and re-emit part of it at lower energies, thus 
performing a spectral down-shift. The radiation is isotropically scattered 
by the luminescent centers and is subjected to refraction upon reaching the 
interface between the slab and the air. Consequently, depending on the 
value of the material's refractive index, part of the radiation is trapped by 
total internal reflection and part of it is lost through the upper and lower 
surfaces. The trapped radiation is wave-guided and reaches the edges of 
the slab, where it can be converted into electricity by standard solar cells. 
These systems do not require tracking and limit the surface of the high-cost 
materials used for the solar cells, thus offering a cheap alternative to stan- 
dard photovoltaic modules. However, to this day, a wide-scale application 
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of these concentrators has been inhibited mainly by their low concentra- 
tion/conversion efficiency, due to several loss mechanisms of which the main 
one is probably the low efficiency of the dyes' wave-guiding properties due to 
low absorption and self-absorption. A comprehensive review on the current 
research status of LSCs is presented in [5]. Different Authors considered 
the problem of modeling these devices, following either a thermodynamic 
approach (based on Chandrasekhar's radiative transfer theory, [^) or a 
computational ray-tracing approach [7]. Both ways represent a broad-scale, 
macroscopic theoretical model of LSCs. In our paper, we focus instead on a 
microscopic model, aimed at characterizing the working principles of the de- 
vice at the molecular scale, considering the interactions between the doping 
dyes and the solar photons. Indeed, the correct physical description under- 
lying the trapping mechanism of a LSC relies on the principles of quantum 
mechanics which govern the molecular spectral conversion process. In doing 
this, we start off from the microscopic thermodynamic model developed by 
Ross [10] and later elaborated by Yablonovitch in his seminal works |1H [T2] . 
With respect to the latter models, we propose a modification of the molecu- 
lar rate equations of the system to correctly account for interaction processes 
involving multiple molecular levels. We explicitly derive the equilibrium ra- 
diation density fields for the molecular system interacting with the solar 
photon bath and describe a method for calculating the parameters involved 
in the relaxation process coupled to the radiative interactions. 

2. Definitions and notations 

For a monocromatic light source S, we define the mean photon-flux density 
at frequency v, ip{v), to be 



(2.1) ^H-/..' 

where h is the Planck constant and / the light intensity. The mean pho- 
ton flux has units photons / cm? x s and may be expressed in terms of the 
radiation energy density per unit volume and unit bandwidth p as 

(2.2) (^(i.) = 

where c has to be substituted by c/n for a material of refractive index n. 
In our discussion we shall be focusing on the interactions between light and 
matter, in the form of molecular or atomic transitions caused by a thermal 
radiation field impinging on a set of molecules. The transitions between 
different molecular energy levels due to electron excitations are related to 
photon absorption and emission. A fundamental physical quantity of our 
model is represented by the transition rate, which expresses the probability 
density per unit time of a transition between atomic energy eigenstates and 
is defined as 

(2.3) W = ip{u)a{u), 
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where W is the transition rate and o"(z^) is the transition cross-section, i.e. 
the effective cross-sectional atomic area (in cm?). As we shaU see later, a 
is determined by the strength of the dipole oscillations which, coupled to 
the incident photon momentum, give rise to the transition. In a LSC the 
light beam inducing molecular transitions comes from the sun and can be 
approximated by a thermal radiation distribution from a black-body source 
of temperature T ~ 5800-fC. The spectral energy density per unit volume 
and unit bandwidth is given by Planck's formula for black-body radiation 

, , . Sirhi/'^ 1 

(2-4) Pi-) = ,3 ,H./KT_y 

The same quantity can be expressed in terms of the wavelength A as 

Svr/ic 1 

(2.5) piX) 



X5 ghc/\KT _i 



3. Interaction dynamics 

As mentioned earlier, whereas the coarse features of a LSC can be treated 
with the principles of geometric optics, the correct dynamics of the molecular 
interactions is described by the principles of quantum mechanics. The model 
developed here is of general validity and can be applied to all cases in which 
fluorescent molecules act as individual spectral converters. In the notation 
of second quantization, the electromagnetic field is an ensemble of photons 
whose state is defined by a momentum and polarization vector for each 
frequency v. The field is specified by giving the number of photons in a 
given state; we set nt^jy to be the number of photons with momentum k and 
frequency v. 

In a volume of space V containing atoms or molecules the energy operator 
(total energy), or Hamiltonian 7i, of the system consists of three terms: the 
radiation field energy Tiem, the molecular energy TLmoi and the interaction 
energy Timt- The latter describes the intensity of the transition processes 
and is used to compute transition rates between different energy levels. We 
introduce the following notation borrowed from Yariv [9]. Let Ek;/ be the 
electric field generated by photons of mode kz/ (polarization and frequency) 
and r the position vector relative to a particle interacting with the field. 
The interaction Hamiltonian can be written as 

'Hint = — eEky • r 

(3.1) 

where e is the electric permittivity, a, are the photon annihilation and 
creation operator of a single radiation mode. In our system, an incident ra- 
diation beam coming from the sun has a distribution of different frequencies 
and momenta, each associated to a specific molecular transition. In general, 
when a photon of mode kz^i is absorbed, the molecular electrons are excited 
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to a higher energy level: the transition energy, which equals the electron's 
levels energy difference, is given by 

(3.2) Ei-Eo = hui. 

We refer hereafter to a general chemical compound that can be schema- 
tized in terms of a three-level system: a fundamental state and two excited 
states relative to absorption and emission. We consider the case in which 
these levels may be treated as well-defined electronic eigenstates. Further- 
more, we make the assumption- common in dealing with fluorescence- that 
spontaneous emission dominates over stimulated emission. 

During absorption, the number of photons in the radiation field decreases 
by one unit (we consider only single-photon processes) while the molecular 
electronic state is shifted to an upper level. If Eq,Ei are the lower and 
higher electronic energy eigenvalues respectively, the rate of absorption is 
given by Fermi's golden rule 

(3.3) Wabs = h\{Eo,ni,^\nint\Ei,n^^ - 6{Ei - Eq - hui), 

where 6 denotes Dirac's delta function. The number of photons in a given 
mode is related to the energy density by 

(3.4) Tik^ = , 

and the average absorption rate is 

27r^e^z/ i2 / \ 

(3-5) Wabs = g^2g l/^l.ol P^^)' 

where /ii^o is the dipole moment of the transition, and ui the associated 
frequency. The same formula expresses also the rate of induced emission, 
where the electronic states return to a lower energy level by emitting one 
photon to the external field. The coefficients of W can be related to Ein- 
stein's coefficients of laser theory as follows. According to Einstein [T3], the 
rate of absorption in the presence of thermal radiation is 

(3.6) W = Bp{ij), 

with p(z/) given by (j2.5p . The constant B is expressed in terms of the above 
coefficients by 

(3.7) ^-l^l""!'' 

and is termed absorption coefficient. For an exhaustive discussion of Ein- 
stein's derivation and its implications on the quantum theory of light-matter 
interactions, see for example Dirac 
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4. Radiative processes 

The theoretical model developed hereafter applies in general to spectral 
conversion procedures in which absorption and emission involve different 
energies and different atomic configurations. In most fluorescent species ab- 
sorption and emission take place in two different parts of the molecule and 
involve different energy levels. The first rigorous thermodynamic description 
of such spectral converters was formulated by Ross [10] and later extended 
by Yablonovitch |1H I12j. However, their model was based on the assump- 
tion that both absorption and emission are fully reversible processes and can 
share chemical equilibrium with the radiation field. As we shall see later, a 
more correct description of a general molecular spectral converter, involves 
at least three, rather than only two, energy levels and therefore is charac- 
terized by a double-equilibrium condition between radiation-absorption and 
radiation-emission. These two processes, both reversible and thus entropy- 
preserving, are related to each other by an internal energy (or charge) trans- 
fer which requires the introduction of two different chemical potentials. As- 
sociated to these processes there is an increase of entropy of the electronic 
system coupled with the photon bath: the entropy acquired upon absorbing 
a photon of energy hi^i is only partially rendered upon emission, as the emit- 
ted photon has an energy hi'2, with f2 < vi- This entropy difference can 
be more rigorously described referring to the change of molecular electronic 
states between absorption and emission by considering the transition from a 
pure electronic state to a statistical mixture of states due to a superposition 
of different nuclear wavefunctions. We can then apply the formula for the 
von Neumann entropy 5 of a quantum system [T7] described by a density 
operator S 

(4.1) S(S) = -TrS log(S), 

where T, is the electronic density matrix corresponding to the emission ex- 
cited states and the log function refers to its diagonal form. Note that this 
analysis assumes that the state of the system before emission is described 
by a superposition of nuclear and electronic wavefunctions of the form 

(4.2) ^' = ^Q,™H,(g)$z,^, 

l,m 

where H, $ denote the electronic, nuclear wavefunctions respectively. With 
the above notation, we have 

(4.3) J:^, = Y,{'^l,m,'f^^^l',m), 

m 

where ^5*^ is the projector on the state (j4.2p . The dynamics of the three 
processes (absorption, emission, thermalization) can be described by a set 
of differential equations for the electronic occupations of the absorption and 
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emission excited states respectively, here labelled as Ni,N2 
(4.4) ^iVi = Bipiui)No - {Bip{ui) + Ai)Ni - qNi 

j^N2 = B2P{u2)Nq - {B2P{V2) + A2)N2 + qNi, 

where B is given by (j3.7p . ^ = and q is the energy transfer intensity, 
specified below. 



Eo 



Figure 1. Sketch of the process dynamics: absorption (be- 
tween levels Eq and Ei), spontaneous emission (between lev- 
els E2 and Eq) and internal energy transfer (between levels 
El and £'2)- Processes occurring with low probability are 
represented by dashed lines and the energy levels refer to the 
electronic states exclusively. 

A is Einstein's coefficient of spontaneous emission, introduced in order 
to account for thermal equilibrium and is simply given by the inverse of 
the spontaneous life-time of the atomic level. The system (|4.5p can be 
solved by considering that, as thermodynamic equilibrium is reached, the 
population fractions referred to the ground state are the relative Boltzmann 
distributions 



(4.5) 



— = exp[-(/ll/i - IJ.l)/KTrr 
N2 

No 



exp[-{hu2 - P2)/kT„ 



where Tm is referred to the molecular temperature after the interaction and 
pi,P2 are the equilibrium chemical potentials relative to z^i,z^2 respectively. 
The set of equations can be solved for the energy densities yielding 

(46) ^,^^,_87rhiyf/c'+q/Bi 



1 
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In the above formulas, we have introduced a frequency-dependent chemi- 
cal potential of the molecular system in equilibrium with the radiation field 
[15j . With regard to the latter quantity (borrowed from the theory of lasers) , 
it is important to observe that from a statistical point of view the molecular 
system is described by a Grand Canonical ensemble in thermal and chemical 
equilibrium with the photon bath of the radiation field. In this case, the 
standard concept of particle exchange is substituted by the more general 
concept of energy quanta exchange. Indeed, the chemical potential of a sys- 
tem constitutes a measure of its ability to exchange particles with a bath at 
assigned values of the other thermodynamic parameters. 

Prom (j4.6p one can note that the presence of the dissipative term q mod- 
ifies the standard energy distributions by increasing the density of the in- 
cident radiation field, while decreasing the emitted one. Thus, the term q 
denotes the rate of energy transfer from radiative to thermal modes due to 
interaction between light and particles. 

This term can be related to the quantum-mechanical description of molec- 
ular dynamics as follows. Our molecular system can be decomposed in three 
constituents: the electronic states involved with absorption, those involved 
with emission, and the remaining nuclear degrees of freedom related to equi- 
librium rotational and vibrational motions. The above electronic states are 
differently localized in the compound; emission states are generally related 
to a metallic ion. 

The total molecular wave-function can be computed within the frame of 
the so called Born-Oppenheimer approximation. The underlying idea of this 
approximation method consists in computing separately the electronic and 
the nuclear wavefunctions, the separation being justified by their different 
kinetic energies. In the first step, the electronic wavefunctions are found by 
solving Schrodinger's equation for fixed values of the nuclear coordinates 



where Ef = Ef{R) are the electronic eigenvalues, dependent on the nuclear 
coordinates R. By varying adiabatically the value of R one can determine 
the potential energy surface Ef{R). In the second step, the nuclear energies 
and wavefunctions are derived in correspondence of a given electronic eigen- 
value by introducing in the Schrodinger equation the kinetic term which was 
initially neglected 



where contains partial derivatives with respect to i?, which is now a 
variable and no longer a parameter. 

Thus, the state of the electrons depends parametrically on the nuclei posi- 
tions, whereas the one of the nuclei is labeled by the electronic eigenvalues. 
On the basis of these considerations, we can write the total approximate 
state of our molecule as a factorized wavefunction 



(4.7) 



ne{r,R)^i{r,R) = E'i^i{r,R) 




[r„ + Ef{R)]<^i^M = El^iAR) 



(4.9) 
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where r, R denote the electronic and nuclear coordinates respectively, H the 
electronic component and $ the nuclear one. 

Upon absorbing a photon of frequency i^, the electronic energy levels of 
the molecule undergo a transition from a state to a state H/^ , with hv = 
Ef^ — Ef^. In parallel, interaction with the incident photons causes a change 
in the nuclear momentum with a relative increase of the kinetic energy, from 
a value Emi to a value Em2- This transition corresponds to a vibrational 
motion of the nuclei away from their original equilibrium positions and a 
rearrangement of the global molecular coordinates. In turn, nuclear motion 
causes a perturbation of the electronic wavefunctions which determines a 
change in charge distribution. The latter may be described as an electron 
transfer which involves a (possibly) non-radiative transition from the original 
level Eij^, acquired after absorption, to an intermediate state Ei^ + 5Ei^, 
before finally reaching the emission state Ei^. The transition rate of this 
internal energy transfer can be estimated from the principles of perturbation 
theory applied to the action of the nuclear kinetic energy operator on the 
unperturbed wavefunctions (j4.9p 

(4.10) q=\^^,Tn^f\^S{E2-Ei), 

where = H/^ (r, i2)(8)$/^^mi i^) total molecular state after absorption, 
= ^.[^{r^R) ® ^>;2,m2(-^) the one before emission, and the delta function 
5{E2 — El) refers to the total energy balance of the transfer including both 
the electronic and the nuclear shifts, with £^1^2 = 2 '^-^1 2' ^^^^ order of 
magnitude of q can be computed using quantum chemical methods, which 
allow to determine the overlap between the two molecular wavefunctions 
involved in the process. 

After reaching the state the electrons return to their original energy 
level by emitting a photon of lower frequency z^', such that hu' = Ef — Ef . 



5. Incident and emitted radiations 

The incident and emitted energy densities are represented by (14.61) . with 
chemical potentials 

(5.1) = E,(^l-^^, 

where Tq,Ti are the ambient temperature and the incident field black-body 
temperature respectively. In a real setting, the incident radiation causes 
thermal energy exchanges between the particles and the photon field. After 
absorption, the system transfers energy to a lower electronic level and thus 
loses kinetic energy in terms of heat by raising its temperature from To to a 
higher value Tm- (Note that this temperature has been already introduced 
above as the correct equilibrium temperature of the molecules). This process 
corresponds to thermalization of the molecular system due to non radiative 
internal energy transfer (the rate of which is given by the term q above). 
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Tjn can be estimated thanks to the equipartition theorem of statistical me- 
chanics, which states that in a system at equilibrium, the temperature is 
equally distributed between the different forms of internal energy. It follows 
that the molecular temperature Tm can be computed by the internal energy 
average using the relation 

(5.2) Erv 

where k is Boltzmann's constant and Erv is the molecular roto-vibrational 
energy. The latter quantity can also be calculated by quantum chemical 
simulations, as well as experimentally measured. The relative numerical 
values and figures are included in a separate file. 

6. CONCENTRATION AND EFFICIENCY 

From (14. 6p . we see that the energy density of the emitted field at the 
emission frequency is higher than the one of the incident field at the absorp- 
tion frequency. On the basis of the above treatment, we could estimate an 
ideal spectral concentration factor directly related to the conversion process, 
defined by the relation 

(6.1) Cm = -, 

Pi 

which, neglecting stimulated emission, yields 

(6.2) Cm = . ^^.,+^,-h.2-hu,/.T^_ 

STrhuf/c-^ + q/Bi 

Ideal efficiencies of LSCs are however limited by the combination of four 
loss mechanisms: absorption efficiency (rja), quantum efficiency of fluores- 
cence (r?/), self-absorption coefficient {9q) and total internal reflection (r/t). 

To give a concrete example, let us consider a slab of refractive index 
n = 1.5 doped with a compound having r]a = 0.05, r/j = 0.7 and 6q = 
0. Defining the geometric concentration factor C as the ratio between the 
upper, absorbing surface and the edge surface, we obtain the total flux gain 
as 

(6.3) G$ = G7?ar//r/t(l - Bg). 

Above we report a plot of the incident solar spectrum as a function of 
wavelength and of the emitted fluorescence spectrum modulated by the 
molecular blackbody radiation which propagates toward the edges of the 
slab. 

7. Conclusions 

In this work we presented a microscopic model for the interaction pro- 
cesses taking place in a LSC. These processes involve multiple molecular 
levels, as they are usually coupled to non-radiative, besides than radiative, 
energy transfer mechanisms. We derive the rate equations for the molecular 
system in equilibrium with the solar photon bath and highlight a way of 
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Figure 2. Energy density vs wavelength of incoming (left) 
and wave-guided (right) radiations. 

determining the parameters involved in the non-radiative processes, which 
are a subject of a further investigation currently in progress. Furthermore, 
we provide examples of the type of spectral distributions involved. In order 
to fully understand the efficiency of LSCs it appears fundamental to han- 
dle not only the macroscopic, ray-tracing aspect of the device, but also the 
molecular-scale efficiencies involved. In these respects, it would be of great 
importance to be able to predict on the basis of molecular properties the 
entity of thermal relaxation and fluorescence efficiency. 

Acknowledgements. P.F.S. would like to thank L.C. Andreani for helpful 
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